Localization of Two-Component Bose-Einstein Condensates in Optical Lattices 
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We reveal underlying principles of nonlinear localization of a two-component Bose-Einstein con- 
densate loaded into a one-dimensional optical lattice. Our theory shows that spin-dependent optical 
lattices can be used to manipulate both the type and magnitude of nonlinear interaction between the 
ultracold atomic species and to observe nontrivial two-componentnlocalized states of a condensate 
in both bands and gaps of the matter-wave band-gap structure. 
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Recent experimental studies of Bose-Einstein conden- 
sates (BECs) in optical lattices, including Bloch-band 
spectroscopy , patterned loading of the lattice , and 
nonlinear tunnelling Q , demonstrate unprecedented lev- 
els of manipulation of the coherent matter waves. Tun- 
ability of the periodic lattice potential combined with in- 
trinsic nonlinearity of BEC provides an excellent testbed 
for many concepts of condensed matter physics and non- 
linear physics. In particular, strong parallels have been 
drawn |^ between nonlinear optics of ultracold atoms in 
optical lattices and behaviour of optical waves in peri- 
odic nonlinear media, where instability and localization 
phenomena have been well explored |5|. Akin to optical 
structures with periodic refractive index, such as waveg- 
uide arrays and photonic crystals, optical lattices form 
a band-gap structure for coherent matter waves which 
modify diffraction properties of atomic wavepackets [ij. 
The coexistence of both normal and anomalous diffrac- 
tion regimes was predicted to lead to a nonlinear localiza- 
tion of condensate with both attractive and repulsive in- 
teractions in the gaps of a linear Bloch-wave spectrum . 
Recent experiments on diffraction management of BEC 
in lattices @i P| have paved the way to the first observa- 
tion of a localized atomic wavepacket (bright gap soliton) 
in a BEC with repulsive interactions [lOj . 

Most of the theoretical and experimental studies of co- 
herent matter waves in optical lattices are concerned with 
single-component BECs. However, it is well known that 
condensate mixtures display many novel physical effects 
not found in single-species BECs including domain 
separation of the ground state and stabilization of non- 
linear localized states. Experiments on spin-dependent 
optical lattices |l2l| have shown that a two-component 
BEC composed of two hyperfine states of the same atomic 
species can be effectively and coherently manipulated in 
an optical lattice. Recently, analysis of periodic (Bloch) 
states of a multi-component BEC in a lattice [l^j has 
revealed that the modulational instability of excited peri- 
odic states can potentially lead to the formation of multi- 
component localized states (solitons) 0. This conclu- 
sion agrees with the recent theoretical study of nonlinear 
inter-gap localization of multi-component coherent light 
fields in periodic optical structures jl5j- 

In this Letter, we address the problem of nonlinear 



localization of a two-component condensate and forma- 
tion of two-component matter- wave solitons in an optical 
lattice. We show that an optical lattice can be used to 
modify the type of the effective nonlinear interactions 
both within and between the condensate species, without 
the Feshhach resonance manipulation of the scattering 
lengths. As a consequence, novel types of nonlinear lo- 
calization of coherent matter waves can be achieved both 
in gaps and bands of the linear Bloch-wave spectrum. 
In particular, when one of the condensate components is 
in the periodic Bloch state, the second component ex- 
hibits effective periodic potential combined with the op- 
tical potential and that induced by the mean-field of the 
Bloch wave. As a result, nontrivial localization of a two- 
component BEC in a self-induced lattice can occur in the 
form of a two-component "dark-bright" band-gap soliton. 
The effects of two-component localization discussed here 
can be explored in the current BEC experiments, as well 
as in the optically-induced photonic lattices |^. 

We model the dynamics of a two-component BEC with 
repulsive interactions in an optical lattice by the mean- 
field equations for the wavefunctions of the condensate 
components \a) and \h): 
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where L„ = -d'^jdx^ + Yn{x) + gnn\^n\ 
gah^a^X- Here we assume that the two 
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^ and Lab = 
BEC compo- 
nents can exhibit different potentials, namely: Vaix) — 
Vb s\T?(khx) and Vb[x) = Vb sin^(A;La; + 6*), which can be 
shifted relative to each other, e.g. by varying the polar- 
ization angle in the case of a spin-dependent lattice 
These potentials are less constrained than those in 
to allow for consideration of a general setup where the 
condensate components may actually belong to different 
atomic species |16ij . being manipulated independently. 
The one-dimensional model is derived for strongly 
anisotropic BEC clouds with the tight confinement direc- 
tion transverse to that of the lattice and made dimen- 
sionless by adopting the lattice units of length xq = k'j} 
(typically ~ 10~^ /zm) and energy = h?k\l1m. The 
stationary states of the wavefunction of the n-th BEC 
component are found as 5'„(a;,i) = ijjn{x) exp(—ifj,„t), 
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FIG. 1: Chemical potential and momentum of the BEC com- 
ponents, relative to the lattice band-gap structure, necessary 
to achieve (a) attractive-attractive, (b) repulsive-repulsive, 
and (c) repulsive-attractive effective nonlinear interactions. 
For the fixed parameters of the \a) component either of the 
|fe) locations can be chosen. Colors code different signs of 
effective diffraction. 



where fin is the chemical potential normalized by Er- 
Without the repulsive mean-field interactions, the ground 
and excited states of each of the BEC components are 
periodic matter waves: ipn = (t>nix)Gxp(ikx), where the 
lattice momentum k lies within the Brillouin zone. The 
periodicity of the lattice leads to the band-gap struc- 
ture of the spectrum At the j-th edge of the 
Brillouin zone /i„ = fin\ the wavefunction (pn^x) = 
(f'ri^x + 7r/fci) is a Bloch state with the lattice peri- 
odicity. Schematics of the Bloch-wave spectrum in the 
reduced-zone representation (|fc| < 1) is shown in Fig. ^ 
A stationary BEC wavepacket located, in the momen- 
tum fc— space, near the j-th band edge can be described 
as a corresponding Bloch wave, dressed by a slowly vary- 
ing (in x) envelope, '^riix,t] fii^'') = ^n{x,t)(l)i^\x). 
Considering that the condensate components \a) and \b) 
can be positioned at the band edges i and j, respectively, 
the envelope equations take the form [isj : 
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where the effective diffraction of the wavepacket is de- 
fined by the curvature of the corresponding band edge, 
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and the effective inter- 



action coefficients depend on the symmetries of the Bloch 
states at the band edges: gj^l^^ = gi^^ J 14)^^ (^l^^ I'^dx. 
The relative shift of the spin-dependent lattice potentials 
can enhance or dampen the cross-species nonlinear in- 
teraction without any change in the effective diffraction, 
which introduces a new degree of experimental control 
into the system. Eqs. Q and their complex conjugate 
form are equivalent, so even for the case of anomalous 
diffraction (_D„ < 0) 'i'l the equations can always be re- 
written in the more intuitive normal diffraction form, but 
with opposite signs of the effective nonlinear coefficients. 

Envelope equations |2Jl predict three types of the ef- 
fective nonlinear interactions that can be realized for 



two-component BECs in a lattice. These regimes can be 
achieved by appropriate placing of the BEC wavepackets 
parameters relative to the band structure, as illustrated 
schematically in Figs.nja-c). Below we consider only the 
main band {Bl: fi^^^\k=o < fJ- < A*^^'^|fe=i) and the first 
excited band {B2: ^^^'^'■|fc=i < fi < fi = A^^''^|fc=o) of the 
lattice that can be populated in experiments []|. 

Attractive- attractive interaction regime [Fig.^a)] can 
be realized when gaa,bb^ 9ab,ba < 0, \gaa\ = \gbb\ = \gab\ = 
\gba\- In this regime, near the band edge fi = f-i^'^^k=i 
both condensate components exhibit anomalous diffrac- 
tion, and therefore can exhibit self-focusing in the form of 
two- component bright gap solitons. The envelope equa- 
tions 101 take the form of the integrable Manakov sys- 
tem, well studied in the context of nonlinear optics (see, 
e.g., According to Eq. (j^J, both localized com- 

ponents can be treated as ground states of the same 
effective potential; they have equal chemical potentials 
Ma = Mb a-nd widths of the spatial density distribution 
envelopes: ^a{x) = $o(a;)cosQ;, ^b{x) = $o(a^)sina, 
where ^q{x) = V Ma/5aaSech(-y2/ia2;), and a is an 
arbitrary parameter. The existence domain for two- 
component solutions in the parameter space {fia, fJ-b} can 
also be calculated using linear waveguiding principles. 
Assuming that the densities of the two component are 
largely dissimilar, one can decouple Eqs. and find, 
numerically, spatial profiles of single-component bright 
gap solitons ipaix] jia) in the entire gap jJ^^ < fJ-a ^ fJ-i^^ ■ 
The low-density second component, ipb can then be found 
as a ground state (i.e. lowest-order guided mode) of 
the effective linear potential combined of the lattice po- 
tential and the mean-field potential induced by the \a)- 
component. The cut-off values i^bilJ-a) for this ground 
state define the lower boundary of the existence domain 
for two-component localized states of the condensate, and 
the upper boundary is /if, ~ fia- The linear waveguiding 
analysis of the complete model Q confirms the predic- 
tion of the envelope theory that in the degenerate case 
of equal effective interaction coefficients the existence do- 
main is confined to the line /ib = fJ-a (see Fig. 2). Re- 
markably, this degeneracy can be lifted by shifting the 
lattices trapping the two condensate components relative 
to each other. The shift destroys the equality between the 
interaction coefficients so that g'ab/gaa > 1 and, accord- 
ing to the envelope theory, yields the new cut-off values 
for theV'& modes, ^6 = (Ma/4)(l-\/l + ^ah/g aa) > l^a- 
The boundaries of the significantly expanded existence 
domain calculated from the guided-mode analysis of the 
complete model (QJ are shown in Fig. [3 by dashed lines; 
such a dramatic expansion can be achieved by the rela- 
tively small shift of = it/ A. 

Repulsive-repulsive regime [Fig. H^b)] can be achieved 
when gaa,bb,gabM > 0, \gab\ = \gba\, and the wavepack- 
ets are located either at the same, \gaa\ — \gbb\, or dif- 
ferent, \gaa\ 7^ \gbb\, band edges. Let us first consider 
the BEC wavepacket located at the edge of the sec- 
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FIG. 2: Existence domain for the bright-bright atomic gap 
solitons in the {^a,^^b} plane {Vo = 2, Qaa = Qbh = Qab = 1). 
Shaded - two lowest Bloch bands Bla,b and B2a,b exhibited 
by the condensate components \a) and |fe). Solid line - ex- 
istence domain = fib predicted by the envelope theory 
and linear waveguiding calculations (see text). Dashed lines 
- borders of the existence domain for the relative lattice shift 
9 = 7r/4. Top panels: examples of the condensate wavefunc- 
tions at the marked points of the existence domain. 



ond band {B2a), Ha = Ma |/c=i- It experiences normal 
diffraction and can support a dark soliton with the en- 
velope $a = \/ Ma/5aatanh(^/iaaa;), imprinted onto the 

extended background of a Bloch wave ^i^^ [x) 0| . Orig- 
inating at the band edge, such dark states with zero 
group velocity can exist within the entire band. Their 
spatial structure suggests that they can be created ex- 
perimentally by phase-imprinting technique, which docs 
not require access to the gap. It is easy to see that, in a 
two-component BEG, the Bloch state in |a)-component 
induces a periodic potential for the |6) component which 
acts together with the potential of the optical lattice. 
As a result, the band-gap structure of the Bloch- wave 
spectrum for the 1 6) -component is significantly modified. 
The original bands are shifted, so that the edge of the 
modified second band /i^ coincides with the chemical 
potential of the nonlinear Bloch state (j)^^ (see Fig. 
The remarkable effect of this condensate-induced lattice 
potential is that the |&)-component can now be spatially 
localized in every gap of the induced band-gap structure 
in the form of a bright gap soliton. Moreover, the ex- 
istence domains for the coupled states calculated using 
the waveguiding properties of the dark in-band solitons, 
lie entirely within the original bands B\ and B2. This 
type of localization could have a dramatic experimental 
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FIG. 3: Existence domain for the dark-bright atomic soli- 
tons (Vb — 2, Qaa = gbb = gab = !)• Shaded - two lowest 
Bloch bands Bla,b and B2a,b exhibited by the condensate 
components \a) and \b). Solid lines - borders of the existence 
domains in repulsive-repulsive regime derived from the linear 
waveguiding calculations (see text). Dotted lines - borders 
of the existence domains for the lattices shift 8 = 7r/4. Top 
(bottom) panels: examples of the condensate wavefunctions 
in the repulsive-repulsive (attractive-repulsive) regime at the 
marked points of the existence domains. 



manifestation, whereby formation of a bright soliton in 
one of the condensate components can be achieved by 
phase-imprinting onto the Bloch state of the complemen- 
tary component in the spectral band. This effect can 
be observed when both components are located either 
at the same or the opposite band edges [Fig. 1(b)], but 
always in the normal diffraction regime. The localiza- 
tion of a single-component repulsive BEG in the form 
of a bright soliton is impossible in this regime. There- 
fore, although the dark-bright state can be dynamically 
stable [see Figs. 0] for a dark state imposed onto both 
modulationally stable [Fig^(c)] and unstable [Fig0|(a)] 
Bloch-wave backgrounds, decoupling of the condensate 
components leads to rapid spread of the bright localized 
state, proportional to its mean-field energy [Figs.^b,d)]. 

Dark-bright solitons supported by effectively repulsive 
nonlinear interactions have been previously discussed in 
the context of nonlinear optics [l8j and BEG in a har- 
monic potential [l9| . The striking difference between 
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FIG. 4: Temporal evolution of a dark-bright localized state in 
the repulsive-repulsive regime with both components within 
(a,b) the second band, at fia = 2.2, Hb = 2.1, and (c,d) the 
first band at fia = 0.81, /it — 0.8. Shown are (a,c) initial 
spatial density profiles and profiles of the coupled state and 
a spreading bright component decoupled from the dark one 
at (b) t = 1.35 mc and (d) t = 2.7 mc. Note the slower 
expansion for the lower density state (c). 



the dark-bright locahzed states in the optical lattice and 
dark-bright BEC solitons in a harmonic potential is that 
the bright component in the lattice is localized in the 
gaps of the induced band-gap structure. Therefore the 
multiple existence domains for the two-component state 
always lie between the cut-off line of the ground mode 
guided by a dark state and the corresponding edge of the 
induced band. In contrast with the case of bright-bright 
gap solitons, these existence domains shrink (dotted lines 
in Fig. 3) with the relative shift between the lattices 
Va{x) and Vb{x), due to reduced effect of the induced pe- 
riodic potential on the bands shift. Another drastic fea- 
ture of the dark-bright states in the lattice is the energy 
difference between the states centered at the lattice min- 
ima and maxima, which inhibits mobility of the localized 
states across the lattice. As a result, the interaction of 
two dark-bright states with zero relative velocity (which, 
without the lattice, dep ends of the relative phase of the 
bright components jl^) is suppressed, and their spatial 
separation can be constant with time. 

Repulsive- attractive regime [Fig. ^c)] can be accessed 
when gaa,gab < 0, > 0, I 9ab\ — l^boli Slid 

\9aa\ 7^ \gbb\- The BEC components should be located at 
different band edges, which is more difficult to arrange 
experimentally. However, in the absence of optical lat- 
tice, this regime could only be explored for a mixture of 
BECs with opposite signs of the scattering lengths. In 
this regime, the envelope equations (01 predict the ex- 
istence of "normal" and "reverse" pairs of dark-bright 
solitons . In the "normal" dark soliton on the 

Bloch-wave background in the normal diffraction regime 
(e.g., at > fJ-a^^) is coupled to a bright gap state in 
the anomalous diffraction regime (e.g., at fii, > ij,^2}b)- 
Both components can exist independently as dynamically 
stable localized states, and the coupling is achieved for 
9ab < \9aagbb\- Typical examples of condcnsatc wavefunc- 
tions in this regime are shown in Fig. 3(c,d). The bright 
localized component \a) exists in the entire gap, whereas 
the dark \b) component forms a coupled state only in the 



vicinity of the band edge. Although in the homogeneous 
case the two-component state is highly unstable with re- 
spect to the mutual displacement of its constituents [lOl , 
the optical lattice potential has a stabilizing effect on its 
dynamics. The "reverse" combination, when the dark 
state in the anomalous diffraction regime (effectively at- 
tractive condensate) is coupled to a bright state in the 
normal diffraction regime (effectively repulsive conden- 



sate) requires that ^^j, > 4|^aa56f)| [201; which is hard to 
achieve for realistic experimental conditions. 

In conclusion, we have analyzed the two-component 
BEC in optical lattices and shown that spin-dependent 
lattices can be used to vary both the type and magnitude 
of nonlinear interactions without manipulating atomic 
scattering lengths. Using a simple envelope theory com- 
bined with the principles of linear waveguiding, we have 
demonstrated that localization of BEC in the form of 
gap solitons can occur both in the gaps and bands of 
the linear Bloch-wave spectrum. The formation of the 
dark-bright composite state is a striking effect that does 
not require accessing the spectral gaps, which makes it 
attractive for possible experimental observations. 

We thank M. Oberthaler for many fruitful discussions. 
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